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Abstract: We investigate the discretized version of the thermodynamic Bethe ansatz 

^ I equation for a variety of 1+1 dimensional quantum field theories. By computing Lyapunov 

\^ ' exponents we establish that many systems of this type exhibit chaotic behaviour, in the 

^^ ■ sense that their orbits through fixed points are extremely sensitive with regard to the 

o 
o 

I ' 1. Introduction 

r~| ! The thermodynamic Bethe ansatz (TBA) equation [|^, Q is an important tool in the context 

of 1+1 dimensional integrable quantum field theories. It serves to extract various types of 
'kj^ \ informations, such as the Virasoro central charge of the underlying ultraviolet conformal 

\-{ ' field theory Q, vacuum expectation values |l|, ^ etc. As it is a nonlinear integral equation, 

it can be solved analytically only in very few circumstances. In general, one relies on 
numerical solutions of its discretised version 

CO 

e;^+i(0) = rmA cosh^-^ j de'^As^O - e')\ii(l + e-'B^.e')'^ . (i,i) 

Here r is the inverse temperature, tua is the mass of a particle of type A, is the rapidity 
and ^AB denotes the logarithmic derivative of the scattering matrix between the particles 
of type A and B. The unknown quantities in these equations are the pseudo-energies 
£a{&)- The standard solution procedure for ( |1.1| ) consists of a consecutive iteration of the 
equation with initial values e^(^) = rTn^iCosh^. At the heart of this procedure lie the 
assumptions that the exact solution is reached for n — > oo, i.e. the sequence converges, and 
furthermore that the final answer is non-sensitive with regard to the initial values, that 
is its uniqueness. In general these assumptions are poorly justified and only few rigorous 
investigations for some simple models exist p, ^, ^. So far the outcome has always been 
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that these assumptions on the existence and uniqueness of the solution indeed hold, albeit 
for certain systems convergence problems in certain regimes have been noted [0, |8|. The 
main purpose of this note is to show that this believe has to be challenged and is in fact 
unjustified for certain well defined theories. We note that our findings do neither effect the 
principles of the TBA itself nor the consistency of the quantum field theory it is meant 
to investigate. However, they indicate that one needs to be very cautious when using the 
above solution procedure and making deductions about the physics for such theories as one 
might just be mislead by the non-convergence of the mathematical procedure used to solve 
the TBA-equations. 

Here we will not analyze the full TBA-equations ( |1.1[) , but rather concentrate on 
the ultraviolet regime, that is r ~ 0, in which it possess some approximation. Clearly, 
the occurrence of chaotic behaviour in this regime will have consequences for the finite 
temperature regime. We encounter the interesting phenomenon that the iterative procedure 
is convergent beyond the ultraviolet (or for a certain choice of parameters in some theories), 
but that unstable fixed points are present in the ultraviolet, meaning that this regime can 
never be reached by the iterative solution procedure for (|1 . ![) . 



2. Stability of fixed points and Lyapunov exponents 

For completeness, let us first briefly recall a few well known basic facts concerning the 
nature of fixed points which may be found in standard text books on dynamical systems, 
see e.g. Q. The objects of our investigations are difference equations of the type 

Xn+l=F{Xn), (2.1) 

where n € Nq, x„ G M and F : M — > M is a vector function. We are especially interested 
in the fixed points Xf of this system being defined as 

F{xf) = Xf . (2.2) 



The fixed point is reached by iterating ( p.l[) , if for a perturbation of it, defined as yn = 
Xn — Xf, we have lim„^oo Vn = 0. From ( |2.1| ) we find 

yn+i + Xf = F{yn + Xf) = F{xf) + J ■ yn + 0(|yn|^) for yk ^ (2.3) 

where J is the i x £ Jacobian matrix of the vector function F{x) 

dFi 



'^ dxj 



forl<i,j<£. (2.4) 



The linearized system which arises from (2^) for yn ^ 



Vn+l = J -Vn (2.5) 

governs the nature of the fixed point under certain conditions [y. Evidently, it is solved 

by 

Vn = QiVi with J ■Vi = QiVi for 1 <i <£ . (2.6) 
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Excluding the case when the eigenvectors Vi of the Jacobian matrix are not linearly inde- 
pendent, we can expand the initial value uniquely 



yo 



E 

i=l 



^iVi 



such that 



Vn 



i=l 



(2.7) 



(2.8) 



It is now obvious from (2^) that the perturbation of the fixed point y„ will grow for 
increasing n if |(7j| > 1 for some i S {1, . . . ,£}. In that case the fixed point Xf is said to 
be linearly unstable. On the other hand, if \qi\ < 1 for all i £ {1, ■ ■ ■ ,^} the perturbation 
will tend to zero for increasing n and the fixed point Xf is said to be linearly stable. It can 
be shown that under some conditions Q the fixed points are nonlinearly stable when they 
are linearly stable. 

In general, that is for any point x rather than just the fixed points xj, stability prop- 
erties are easily encoded in the Lyapunov exponents Aj. Roughly speaking the Lyapunov 
exponents are a measure for the exponential separation of neighbouring orbits. One speaks 
of unstable (chaotic) orbits if Aj > for some i £ {!,...,£} and stable orbits if Aj < 
for all i G {!,...,£}. For an arbitrary point x the i Lyapunov exponents for the above 
mentioned system ( p.l[) are defined as 



A,; 



lim 

n— >oo 



n 



Q^[F''m 



lim 

ra— >oo 



^ n— 1 

-Vln 
n ^^ 



Q^[F\x)] 



k=0 



(2.9) 



where the qi (x) are the eigenvalues of the Jacobian matrix as defined in ( p.6D , but now at 
some arbitrary point x. Taking the point to be a fixed point, we can relate ( |2.9| ) to the 
above statements. At the fixed point we have of course F^{xf) = Xf, such that 



A,; 



lim 

n— >oo 



^ n—1 

- Vln|gi(f/)| 
n ^-^ 



k=o 



Inlo'j 



(2.10) 



Therefore, a stable fixed point is characterized by Aj < or |gj| < 1 for alH G {1, . . . ,^} 
and an unstable fixed point by Aj > or |gj| > 1 for some i G {!,...,£}. We can now 
employ this criterion for some concrete systems. 

3. Unstable fixed points in constant TBA equations 



We adopt here the notation of [^, 11, |12[, by which a large class of integrable quantum 
field theories can be referred to in a general Lie algebraic form as g|g-theories. Their 
underlying ultraviolet conformal field theories can be described by the theories investigated 
14, |l5| (and special cases thereof) with Virasoro central charge c = iih/{h+h). Here 



in Ic 



i{i) and h{h) are the rank and the Coxeter number of g(g), respectively. In particular, 
g| Ai is identical to the minimal affine Toda theories (ATFT) jl^, 17| and A„|g corresponds 
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to the g„+i -homogeneous Sine-Gordon (HSG) models [18, |l^. In this formulation each 
particle is labelled by two quantum numbers (a, i), which take their values in 1 < a < £ and 
1 < i < i. Hence, in total we have i x i different particle types. It is a standard procedure 
in this context ||l], |2| to approximate the pseudo-energies in ( [1.1| ) by £^{0) = e\ = const 
in a large region for 9 when r is small. For convenience one then introduces further the 
quantity x\ = exp(— e^) such that ( |1.1| ) can be cast into the compact form 

i I 
< = n 11(1 + x^,)Ki =: Fi,{x) with iV:^ = 5a,5,, - R-^'K,, . (3.1) 

6=1 j=l 

The matrix A'^'^ in (^]^) encodes the information on the asymptotic behaviour of the scat- 
tering matrix. As stated in ( |3.1| ) it is specific to each of the g|g-theories with K and K 
being the Cartan matrix of g and g, respectively. The equations (3.1) are referred to as the 



constant TBA-equations. They govern the ultraviolet behaviour of the system and their 
solutions yield directly the effective central charge 

a=l 1=1 ^ "' 



with C{x) = X^,^]^ x^ jr? -|- ln2;ln(l — 2;)/2 denoting Rogers dilogarithm (see e.g. |20| for 
properties). 



"^ " dxi 



Let us now discretise (3J) and analyze it with regard to the nature of its fixed points. 
According to the argument of section 2, we have to compute first of all the Jacobian 
matrices of F 

= <iT^%r • ^'-'^ 

Next we need to determine the eigensystem of the Jacobian matrix. This is not possible to 
do in a completely generic way at present, since not even the solutions, i.e. fixed points, 
of ( |3.lD are known in a general fashion. Instead, we present some examples to exhibit the 
possible types of behaviour. 

3.1 Stable fixed points 

We start with a simple example of a stable fixed point. We present the ^2|^i case, which 
after the free Fermion (j4i|Ai) is the next non-trivial example in the series of the minimal 
ATFTs, the scaling three-state Potts model with Virasoro central charge c = 4/5. The 
TEA has been investigated in H, |2[. The constant TBA-equations ( ^.1| ) 

XI = (1 + Xi)-1/3(1 + ^2)-2/3 ^ ^^(^)^ X2 = (1 + Xi)-2/3(l + X2)-^'^ = F^ix), (3.4) 

can be solved analytically by the golden ratio r := (\/5 — l)/2 = xi = X2- Using this 
solution for the fixed point xj = (xi,X2), we compute the Jacobian matrix of -F(x) 

1 / 2 2T-2 \ 
J(*/) = -3 2,. ,. . (3.5) 
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with eigensystem 



qi 



-r^ PS -0.38197, ^i = (l,l), g2 = rVS « 0.12732, t;2 = (-l,l) 



As the eigenvectors are obviously hnearly independent and \qi\ < l,\q2\ < 1, we deduce that 
aU Lyapunov exponents are negative and therefore that the fixed point is stable. Indeed, 
numerical studies of this system exhibit a fast convergence and a non-sensitive behaviour 
with regard to the initial values of the iterative procedure. 

For some theories the solutions are known analytically in a closed form. For instance, 
the constant TBA equations for the j4i|^£-theories (= SU{i + l)2-HSG-niodel) are solved 

by 



X 



sin[7r(l + i)A] 



1, 



for 1 < i < 



(3.6) 



with A = 1/(3 + £) 
Jacobian matrix (| 



sin(7rA) 

23, |l^. Taking this solution for the fixed point we compute the 
with Nii = (5jj+i + (5jj_i)/2 to 



•y 



jiJt^ 



'11 



(%) 



sin[7r(2 + i)A] 
sin(i7rA) 



^i,j+l + 



sin(i7rA) 
sin[7r(2 + i)A] 



^«j- 



(3.7) 



and the eigenvalues to qi = cos[7r(i + 1)A]. As \qi\ < 1 for all i G {1, . . . ,i} the fixed points 
are stable. We investigated various minimal affine Toda field theories which all posses fixed 
points of this nature. In fact, the general assumption is that all systems exhibit such a 
behaviour. We present now some counter examples which refute this believe. 



3.2 Stable two-cycles 

We start with a system which does not possess a stable fixed point, but rather a stable 
two-cycle, i.e. a solution for 

G{x) := F{F{x)) = X . (3.8) 

We consider the A2|^2-theories (= 5'C/(3)3-HSG model) studied already previously by 
means of the TBA in [24|. Its extreme ultraviolet Virasoro central charge is c = 2. The 
constant TBA-equations for this case read 



1 _ (l+xf)^/3(l + xi)V3 
""^ ~ (l+xl)V3(l + ^1)2/3 

2_ {l+x\f/^il + xl)'/^ 
^' ~ (l+x2)l/3(l + ;r2)2/3 

with analytic solution x{ = Xg = 
compute the Jacobian matrix for F{x) 



Ji^f) = 6 



= i"l(^), 


X'2 


= F3{X), 


4 


X 1 — Xo — 


1. ' 



(l + xf)V3(l+ ^2)2/3 
(l + xl)2/3(l+xl)V3 

{I + x\)^/Hl + xlf/^ 

(l + x2)2/3(l+2;2)l/3 



F2ix), (3.9) 

F4(x), (3.10) 



1. Taking this solution as the fixed point, we 



/-I -2 2 l\ 
-2-1 1 2 
2 1-1-2 
V 1 2-2-1 



(3.11) 
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ith e 


gensystem 




qi 


= -1, 


Vl 


Q3 


= 0, 


V3 



(-1,-1,1,1) 

(1,0,0,1), 



92 
94 



1/3, 
0, 



V2 



(-1,1,-1,1), 

(0,1,1,0). 



(3.12) 



We observe that there is one eigenvalue with \qi\ = 1, which is generahy cahed a marginal 
behaviour, i.e. the stability properties depend on the other eigenvalues and on the next 
leading order. In fact, we can see from (|2.8| ) that the perturbation of the fixed point will 
remain the same even for large values of n, flipping between two values and thus suggesting 



the existence of a stable two cycle (|3.8|). We find that (3.8) can be solved by 



x\ 



X2 



l/xl 



1/xl 



(3.13) 



for any arbitrary value of k. To determine the stability of the two-cycle we have to compute 
the Jacobian matrix for G{x) 



J{Xf) 



+ 9k 



/ 5k 4k; -4k^ -5k^ \ 
Ak 5k — 5k^ — 4k^ 
-4/k -5/k 5 4 
\-5/k -4/k 4 5/ 



(3.14) 



which has eigensystem 

qi = l, vi = {-k^,-k'^,1,1) 

93 = 0, V3 = {k, 0,0,1), 



92 
94 



1/9, 
0, 



V2 



-K^,K^,-l,l] 



(0,K,1,0). 



(3.15) 



We conclude from this that one approaches a stable two-cycle when iterating the discretised 
version of (|3!^ ). Thus, we note that the TBA-system for the 5L'"(3)3-HSG model in the 
ultraviolet regime does not posses a stable fixed point but an infinite number of stable 



two-cycles of the type ( 3.15 ). It is now intriguing to note that when using this solution 
to compute the effective Virasoro central charge ( ^.21) , one always obtains the expected 
value c = 2 for any value of k G M, simply due to an identity for the Rogers dilogarithm 
£(1 — x) + C{x) = 7r^/6. Hence, despite the fact, that one is using entirely wrong pseudo- 
energies, one obtains by pure luck an apparent confirmation of the theories consistency. 



3.3 Unstable fixed points, chaotic behaviour 

In this section we present some TBA-systems for well-defined quantum field theories, which 
exhibit a chaotic behaviour in the sense that their iterative solutions are extremely sensitive 
with regard to the initial values. 



3.3.1 Ai\Ai 

This model is the 5'L'^(5)5-HSG model with extreme ultraviolet Virasoro central charge 
c = 8. To reduce the complexity of the model, we exploit already from the very beginning 



the Z2-symmetries in the ^4-Dynkin diagrams and identify x\ = x\ 



Xa 



X|, X2 
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Xn 



X3, X2 



4 
X2 



X3, Xi 



xf 



X4. With these identifications the constant 



TBA-equations can be brought into the form 



2\2 



^^ (1+X{)(1 + X^)2 

, _ {l + x\){l + xl) _ 



= Fi{x)., 
Fsix), 



{l + xi){l + xi 

~~ (l+x})2(l + xl)3 

{l + x\f{l + x^f _ 



--F2{x)., 
F^{x). 



We can solve these equations analyticaUy by one, r and f := l/r = (\/5 + l)/2 



(3.16) 
(3.17) 



™2 
X2 



X4 



X4 
Xi 
X2 



X2 



T. 



X2 



J2 
2^3 



X3 



1, 



(3.18) 



With this solution for Xf at hand we compute the Jacobian matrix of F{x) in ( 3.16| ), (|3.17| ) 



/-1/2 -2r 



J{Xf) 


= 


— r - 
?/2 


-3t2 

1 




V 1/2 


2r 


Now we find the eigensystem 






gi = -T^^- 


-2.6180, 




ga = 4t - 2 ? 


b 0.47214, 




53 = 0, 






g4 = 0, 









1/2 \ 



-f /2 



(3.19) 



-1/2/ 



^71 = (-1,-1, 1,1), 
t;2 = (-l,r2,-f2,l), 
i;3 = (1,0,0,1), 
t;4 = (-2r2,T,l,0). 



(3.20) 



As \qi\ > 1 and the eigenvectors are linearly independent, we deduce that the Lyapunov 
exponent Ai is positive and therefore that the fixed point ( |3.18| ) is unstable. Indeed, 
numerical studies of this system exhibit that any small perturbation away from the solution 



( 3.18 ) will lead to a divergent iterative procedure. 

Nonetheless, by some manipulations of ( |3.16|) , (|3.17| ) one can find equivalent sets of 
equations which posses stable fixed points and can be solved by means of an iterative 
procedure. For example, when simply substituting x\ in F2{x) we obtain the equations 



x\ 



^2 

Xi 



Fi{x) 
F3{x) 



F[{xl 

Fm, 



X2 



X2 



(l + xl)(l + xi) 
Fi{x)=Fi{x), 



F^{S), 



(3.21) 
(3.22) 



which are of this kind. Now all Lyapunov exponents resulting from the Jacobian matrix for 
F'{x) at the fixed point ( 3.1^ ) are negative. One should note, however, that even though 
F'[x) it is easily constructed by trial and error from F{x) for the constant TBA-equations, 
the equivalent manipulations on the full TBA-equations ( |l.lD are quite unnatural, albeit 
not impossible to perform once ( ^.1| ) is analyzed. 

One of the distinguishing features of the HSG-models is that they contain unstable 
particles in their spectrum, whose masses are characterized by some resonance parameters 
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Gij with 1 < i,i < ^. We can now interpret these parameters as bifurcation parameters as 
common in the study of chaotic systems and investigate the nature of the fixed points when 
these parameters are varied. In |^ a precise decouphng rule was provided, which describes 
the behaviour of the theories when some of the cr's become large and tend to infinity. For 
S'C/(5)5 we have for instance the following possibilities 

lim 5C/(5)5 = SC/(2)5 SU{A)^ or lim SU{b)^ = SU{?,)^ ® SC/(3)5 . (3.23) 

(T12— >00 (723— >00 

For the algebras involved we found that the fixed point of SU{A)^ is unstable, whereas the 
fixed points of S'C/(3)5 and SU{2)^ are stable. For our SU^b)^ example this implies that the 
fixed point in SU^S)^ SU{3)5 will be stable, whereas the fixed point in S'[/(2)5 ® S'C/(4)5 
will be unstable. In general, we find that, while approaching the ultraviolet from the 
infrared, once the nature of the fixed point has changed from stable to unstable it remains 
that way. This behaviour can be encoded naturally in standard bifurcation diagrams, which 
we present elsewhere. 

We also found unstable fixed points for other HSG- models related to simply laced 
algebras and g|g-theories which are neither HSG nor minimal ATFT. A priori the behaviour 



is difficult to predict, e.g. whereas -D4I-D4 (see |10|] for the solution of (|3.lD ) and D4\A4 



have unstable fixed points, the fixed point in D4\A2 is stable. 

3.3.2 Ai\C2 

This model is the simplest example of an HSG model related to a non-simply laced algebra, 
namely the S'p(4)2-HSG model with central charge c = 2. In the TEA analysis carried out 
in 0] convergence problems in the ultraviolet regime were already commented upon. In 
fact, we find here that all HSG models which are related to non-simply laced Lie algebras 
posses unstable fixed points. The constant TBA-equations for ^i|C2 read 



X 



1 _ . //I _L ^2^/1 _L ^2^n _L ^2^ - T?. (^\ ^2 _ ^ _ / (^ + ^1) 



} = J(l + xDH + xi)(l + xi) = Fi(f), xi = —-^J ^ ; . ^/ . = F2(x),(3.24) 



"'"(i+x?)(i+xi)(i+xi)-''^^"^^' "'- (rT^V(i+x?K?+xi) -''^('')' ('•'') 

with solutions 

x\ = 3, xl = xl = 2/3 and xl = 4/5. (3.26) 

The corresponding Jacobian matrix for F{x) reads 

/ 9/10 5/3 9/10 \ 

1/12 -1/5 -10/27 -1/5 

1/5 -12/25 -4/9 -12/25 

\l/12 -1/5-10/27 -1/5/ 

with eigensystem 

gi ^ -1.3647, vi = (-3.53, 1, 1.8104, 1), 

g2 ~ 0.3973, V2 = (-80.2656, 1, -20.2124, 1), 

gs « 0.1229, V3 = (2.1956, 1, -0.9180, 1), 

g4 = 0, i;4 = (0,-1,0,1). 



J{Xf) 



(3.27) 



(3.28) 
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Since |gi| > we find a positive Lyapunov exponent and therefore the fixed point Xf in 



( 3.26 ) is unstable. 

We also checked explicitly Ai\Cs,Ai\C4, Ai\G2, A2\B2 and found a similar behaviour. 
Based on these examples we conjecture that the constant TBA-equations ( |3.lD related 
to g-HSG-models with g non-simply laced have unstable fixed points. 

4. Conclusions 

We showed that the discretised TBA-equations for many well-defined quantum field theories 
exhibit chaotic behaviour in the sense that their orbits are extremely sensitive with regard 
to the initial conditions. In particular, we found several examples for HSG- models and g|g- 
theories which are neither HSG nor minimal ATFT. Apart from the statements, that all 
j4i|A£-theories have stable fixed points and apparently all HSG-models which are related to 
non-simply laced models have unstable fixed points, we did not find yet a general pattern 
which characterizes such theories in a more concise way. 

Our findings clearly explain the convergence problems reported upon earlier in [0] and 
we stress here that they do neither effect the consistency of the quantum field theories nor 
the validity of the principles underlying the TBA, but only point out the need to solve these 
theories by alternative means. The closest would be to alter the iterative procedure for 
( |1.1| ) as indicated in section 3.3.1 for the constant TBA-equations, by defining equivalent 
sets of equations which have stable fixed points. Unfortunately, we can not settle with 
these arguments the convergence problems for the models studied in |^], as for those the 
fixed points are situated at infinity. 

Our results clearly indicate that one can only be confident about results obtained from 
iterating ( |1.1D if the nature of the fixed points is clarified. 

Acknowledgments: This work is supported in part by the EU network "EUCLID, Inte- 
grahle models and applications: from strings to condensed matter" , IIPRN-CT-2002-00325. 

References 

[1] A. B. Zamolodchikov, Thermodynamic Bcthe ansatz in relativistic models: Scaling 3-state 
Potts and Lee- Yang models, Nucl. Phys. B342, 695-720 (1990). 

[2] T. R. Klassen and E. Melzer, The Thermodynamics of purely elastic scattering theories and 
conformal perturbation theory, Nucl. Phys. B350, 635-689 (1991). 

[3] A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, Infinite conformal symmetry in 
two-dimensional quantum field theory, Nucl. Phys. B241, 333-380 (1984). 

[4] O. Castro- Alvaredo and A. Fring, On vacuum energies and renormalizability in integrable 
quantum field theories, Nucl. Phys. B687, 303-322 (2004). 

[5] C.-N. Yang and C. P. Yang, Thermodynamics of one-dimensional system of bosons with 
repulsive delta function interaction, J. Math. Phys. 10, 1115-1122 (1969). 

[6] A. Fring, C. Korff, and B. J. Schulz, The ultraviolet behaviour of integrable quantum field 
theories, affine Toda field theory, Nucl. Phys. B549, 579-612 (1999). 



Chaos in the thermodynamic Bethe ansatz 



\7] O. A. Castro- Alvaredo, J. Dreifiig, and A. Fring, Integrable scattering theories with unstable 
particles, The European Physical Journal C35, 393-411 (2004). 

[81 O. Castro- Alvaredo and A. Fring, Constructing infinite particle spectra, Phys. Rev. D64, 
0850051-0850057(2001). 

[91 P. G. Drazin, Nonlinear Systems, CUP, Cambridge, 1992. 

[101 A. Fring and C. Korff, Colour valued scattering matrices, Phys. Lett. B477, 380-386 (2000). 

[Ill C. Korff, Colour valued scattering matrices from non simply-laced Lie algebras, Phys. Lett. 
B501, 289-296 (2001). 

[121 O. A. Castro- Alvaredo and A. Fring, Scaling functions from q-deformed Virasoro characters, 
J. Phys. A35, 609-636 (2002). 

[131 D- Gcpner, New conformal field theories associated with Lie algebras and their partition 
functions, Nucl. Phys. B290, 10 (1987). 

[141 G. V. Dunne, L Halliday, and P. Suranyi, Bosonization of parafermionic conformal field 
theories, Nucl. Phys. B325, 526 (1989). 

[151 J- M. Camino, A. V. Ramallo, and J. M. Sanchez de Santos, Graded parafermions, Nucl. 
Phys. B530, 715-741 (1998). 

[161 A. V. Mikhailov, M. A. Olshanctsky, and A. M. Perelomov, Two dimensional generalized 
Toda lattice, Commun. Math. Phys. 79, 473 (1981). 

[171 D. L Olive and N. Turok, Local conserved densities and zero curvature conditions for Toda 
lattice field theories, Nucl. Phys. B257, 277 (1985). 

[181 Q--H. Park, Deformed coset models from gauged WZW actions, Phys. Lett. B328, 329-336 
(1994). 

[191 ^- ^- Fernandez-Pousa, M. V. Gallas, T. J. HoUowood, and J. L. Miramontes, The 

symmetric space and homogeneous sine-Gordon theories, Nucl. Phys. B484, 609-630 (1997). 

[201 A. N. Kirillov, Dilogarithm identities. Prog. Theor. Phys. Suppl. 118, 61-142 (1995). 

[211 A. Kuniba, Thermodynamics of the Uq{Xr{l)) Bethe ansatz system with q a root of unity, 
Nucl. Phys. B389, 209-246 (1993). 

[221 A. Kuniba and T. Nakanishi, Spectra in conformal field theories from the Rogers 
dilogarithm. Mod. Phys. Lett. A7, 3487-3494 (1992). 

[231 A. Kuniba, T. Nakanishi, and J. Suzuki, Characters in conformal field theories from 
thermodynamic Bethe ansatz, Mod. Phys. Lett. A8, 1649-1660 (1993). 

[241 O. A. Castro-Alvaredo, A. Fring, C. Korff, and J. L. Miramontes, Thermodynamic Bethe 
ansatz of the homogeneous sine-Gordon models, Nucl. Phys. B575, 535-560 (2000). 



10 



